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Abstract: The combination of stable isotope labeling of amino acids in mammals (SILAM) and laser capture microdissection (LCM) 
for selective proteomic analysis of the targeted tissues holds tremendous potential for refined characterization of proteome changes 
within complex tissues such as the brain. The authors have applied this approach to measure changes in relative protein abundance in 
ventral tegmental area (VTA) of the rat brain that correlate to pharmacological perturbations. Enriched "C^'^Nj-lysine was introduced in 
vivo via diet. These animals were sacrificed during the middle of the 12-hour light period to extract isotopically "heavy" proteins, which 
were then used as a reference for extracts from dosed, unlabeled rats. Animals were administered an orexin peptide (Ox-B), an orexin 
receptor antagonist (ORA), or a mixture of both (Ox-B + ORA). All samples were obtained at same phase of the sleep cycle. Labeled- 
pair identification and differential quantitation provided protein identification and expression ratio data. Five proteins were found to 
exhibit decreased relative abundance after administration of an ORA, including a-synuclein and rat myelin basic protein. Conversely, 
six proteins showed increased relative abundance upon antagonist treatment, including 2',3'-cyclic nucleotide 3'-phosphodiesterase. 
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INTRODUCTION 

Biochemical characterization of complex biology and 
diseases is a common challenge facing scientists in the 
biomedical and pharmaceutical sciences. A reductionist 
approach, breaking the system into more basic elements 
and studying the interaction of those elements, is one 
method towards understanding the nature of complex 
biology. In vitro models, although powerful and often- 
times predictive, are at many times insufficient and lack 
an appropriate level of contextual specificity for particular 
areas of biomedical research. For example, neuropathic 
pain and modulation of the sleep/wake cycle are two areas 
that are difficult to investigate based on in vitro models. 
Indeed, these and many other biochemical processes may 
be dynamically entwined with the intracellular, extrace- 
llular, and organ- and organism-level environments and 
require a system-wide approach. 



Increased levels of orexin neuropeptides promote and 
maintain arousal and wakefulness, whereas a reduction in 
orexin signaling promotes sleep (for review, see (Mieda 
& Sakurai, 2013)). The activity of orexin neurons oscil- 
lates over a 24-hour period, with orexin peptide release 
rising during wakefulness, peaking late in the active phase 
and rapidly decreasing throughout the normal inactive 
phase (Gotter et al., 2013; Taheri et al., 2000; Zeitzer 
et al., 2003). Pharmacological antagonism of the orexin 
system affords a novel approach for treating insomnia and 
has demonstrated efficacy in rodents, dogs, and humans, 
including patients with insomnia (Winrow et al., 2011; 
Brisbare-Roch et al., 2007; Hoever et al., 2013; Herring 
et al., 2012; Bettica et al., 2012). The ability to regulate 
biologically important patterns, such as the sleep/wake 
cycle, with predictable periodicity is typically associated 
with the suprachiasmatic nuclei region of the mammalian 
brain. However, recent genetic evidence revealed that the 



Received 28 October 2013; accepted 10 January 2014. 

Address correspondence to John J. Renger, Merck Research Laboratories, 770 Sumneytown Pike, PO Box 4, West Point, PA 19454, 
USA. E-mail: john_renger@merck.com 



136 



LCM and SILAM Proteomic Profiling of Sleep Nuclei 



137 



ventral tegmental area (VTA) also displays activity that 
correlates to the normal sleep cycle in mammals (Winrow 
et al., 2009). Physiologically, the VTA serves as the cen- 
tral region of the mesocorticolimbic dopamine system, 
which regulates cognition, motivational processing, and 
reward impulses. 

Recent advances in mammalian tissue proteomics, 
such as in vivo stable isotope labeling amino acids in 
cells (SILAC; Kruger et al., 2008), have enabled accurate 
relative quantitation within rodent model systems (a.k.a. 
"SILAC Mouse" and stable isotope labeling amino acids 
in mammals [SILAM]; McClatchy et al., 2007). These 
approaches incorporate selective labeling of one or more 
amino acid residues in vivo. The result is the generation 
of endogenously translated, chemically identical, yet 
molecular weight-modified, labeled proteins in tissue that 
can be mixed in reference to a control or perturbed state 
for differential analysis through quantitative mass spec- 
trometry. Examples of sample preparations that incorpo- 
rate such heavy-isotope labeling schemes from diverse 
tissue types have been reported, including those for neu- 
ronal applications (Ishihama et al., 2005) and experiments 
that introduced '^N-label into Sprague-Dawley rats via 
enriched food pellets (Wu et al., 2004; McClatchy et al., 
2007). Further optimization of these results indicated that 
progressive labeling of second-generation rats yielded as 
much as a 94% efficiency among all tissues examined 
(McClatchy et al., 2007). The use of a more selective 
'^Cj-lysine in second- to fourth-generation mice demon- 
strated even further improved labeling efficiency while gen- 
erating less complex mass spectra (Kruger et al., 2008). In 
addition, high-precision surgical technologies coupled with 
laser capture microdissection (LCM) now make it possible 
to isolate targeted nuclei with single-cell precision from 
surgically removed tissue. This approach has been used to 
successfully characterize the protein content or differential 
genomic profiles of a number of specific cell populations 
and organ substructures (Liotta & Petricoin, 2000; Geho 
et al., 2004; Sanders et al., 2008; Patel et al., 2008). 

Here, we apply "C^'^N^-lysine SILAC labeling in 
combination with targeted tissue selection of the VTA and 
a multiarm dosing method for investigating the specific 
biochemical activity of a recently described orexin recep- 
tor antagonist ("compound 18," referred to here as ORA; 
Mercer et al., 2013) in rat models. This antagonist exhib- 
its ~ 1400-fold binding selectivity and a 24-fold functional 
selectivity for orexin receptors OX2R over OXIR, with 
> 100-fold selectivity in a functional screen of 165 other 
targets, and significantly reduces active wake and promotes 
both rapid eye movement (REM) and non-REM sleep in 
rat electroencephalographic (EEG) studies (Mercer et al., 
2013). Our group and others have demonstrated that struc- 
turally distinct 0X2R antagonists effectively promote 
sleep in rodents and that administration of Ox-B peptide is 
sufficient to promote wakefulness (Dugovic et al., 2009; 



Mercer et al., 2013; Winrow et al., 2010). Animals were 
given either the activating Ox-B peptide to promote wake- 
fulness, the orexin receptor antagonist (ORA) to promote 
sleep or a combined treatment (Ox-B + ORA). Relative 
changes in peptide abundance were measured for labeled- 
pairs and then filtered based on the changes that corre- 
late to locomotor activity, including increases relative to 
control with stimulant, decreases relative to control with 
antagonist, or intermediate to peptide and antagonist with 
combined treatment (or the opposite). 

MATERIALS AND METHODS 

Food Preparation, In Vivo Labeling, and Tissue 
Collection 

Unmodified control rat diet and "C^'^N^-lysine modified 
diet (100 g labeled Lys per 10 kg diet) were prepared 
by and purchased from Harlan-Teklad (Boston, MA), 
shipped at — 70°C, and kept frozen upon arrival. ''C^'^N^- 
lysine (> 98%) used in the modified diet was purchased 
from Sigma-Aldrich (St. Louis, MO). The in vivo labeling 
method was followed as reported by McClatchy et al. (2007 ; 
Protocol 2). Briefly, 2-day pregnant Sprague-Dawley rats 
(Taconic, Germantown, NY) were provided ad libitum 
access to water and 25 g of diet per day. Two dames were 
fed unmodified lysine (light) control diet and two dames 
were fed diet enriched with '^C^'^N^-lysine. The dames 
fed the unmodified diet were designated (light dame) LDl 
and LD2 and the dames fed the heavy-isotope enriched 
diet were designated (heavy dame) HDl and HD2. The 
pregnant dames were subjected to a 12 hours lights on 
(6:00 am) / 12 hours fights off (6:00 pm) schedule up until 
weaning, and the pups were continued on this schedule 
and their respective diets for up to seven additional weeks 
prior to necropsy. All pups were carefully inspected and 
deemed healthy. After birth, dame feed portions were 
increased to 50 g/day for 3 weeks. At weaning, males 
were selected to reduce potential experimental complexity 
due to the female reproductive cycle. Each weanling rat 
was housed individually and fed approximately 15 g/day 
of their respective diet for the next 4 weeks. At the end 
of 7 weeks, rats were euthanized for tissue preparation 
as described below. All animals used in this study were 
handled and housed in agreement with the federal animal 
welfare guidelines. All studies were approved in advance 
by the Institutional Animal Care and Use Committee at 
Merck Research Laboratories at West Point. 

Beam Break Activity Measurement and Compound 
Administration 

Male Sprague-Dawley rats (250-280 g) with surgically 
implanted intracerebroventricular (IC V) cannulas (Taconic 
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Farms, surgical model no. SU009; lateral ventricle) were 
individually housed upon delivery. Rats had free access to 
food and water and were maintained on a 12 hours lights 
on (6:00 am) / 12 hours hghts off (6:00 pm) schedule. 
Approximately 24 hours prior to dosing, each animal was 
placed within an infrared photobeam array that tracked 
activity in a xv-plane (ambulations and fine movements) 
and in the z-plane to record rearing behavior (Columbus 
Instruments, Columbus, OH). Baseline activity was 
recorded after all animals were placed in the array. Activity 
data (total beams interrupted) were analyzed using a 
customized Excel template and calculated as area under 
curve. Rats were placed into one of three treatment groups: 
Ox-B (RSGPPGLQGRAQRLIQASGNHAAGILTM, 
Ala", D-Leu'^ C-terminal amide; Tocris Bioscience, 
Ellisville, MO) (2 mM), ORA (30 mg/kg), or Ox-B + ORA. 
The treatment groups were established as follows: on the 
day of the study, the orexin antagonist or vehicle was 
administered orally starting at -10:30 am and dosing was 
staggered at 7-8-minute intervals. After dosing, rats were 
returned to their home cage in the beam break array. Sixty 
minutes after ORA or vehicle dosing, each rat received 
an ICV infusion of either xxx cerebrospinal fluid artificial 
CSF (aCSF) (CMA Microdialysis, North Chelmsford, 
MA) or aCSF containing Ox-B using a CMA Microdialysis 
pump delivering 5 jlL over 60 seconds. The injector can- 
nula was left in place for an additional 30 seconds to 
allow dispersion into the CSF. Rats were then returned 
to the beam break array where the locomotor activity 
was recorded for another 60 minutes. At the end of the 
60-minute period, rats were euthanized by decapitation, 
brains removed, rinsed in phosphate-buffered saline, and 
rapidly frozen on dry ice. Brains were isolated at 7 hours 
after "lights on" ± 30 minutes, representing the middle of 
the light period, and stored at — 80°C for future studies. 



Laser Capture Microdissection and Protein 
Extraction 



The process for isolation of specific brain regions is 
shown schematically in Figure la-e. Brain tissue sec- 
tions (20 jim diameter) were prepared on polyethylene 
napthalate (PEN)-membrane slides (Molecular Devices, 
Sunnyvale, CA) from SILAC-labeled or compound- 
treated rats as described above. Specific brain regions 
were identified by cressyl violet staining using the LCM 
Staining Kit (Molecular Devices) (Figure Ic). Figure Id 
shows the visualized regions that were then targeted 
and isolated by automated laser microdissection using 
the VERITAS platform (Molecular Devices), optimized 
for single-cell accuracy (Rodriguez et al., 2008). Two 
VTA sections were mounted on an individual CapSure 
Macro LCM Cap (MacroCap; Molecular Devices) for 
an average volume of 0.13 mm^ of VTA tissue per cap 
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Figure I. Isolation of specific brain regions from labeled and 
unlabeled rats, (a) Male Fl generation rats fed either control 
or isotopically labeled lysine-enriched chow (see Materials 
and Methods), (b) Rat brain regions are surgically isolated, 
(c) Ventral tegmental area (VTA; outlined in yellow) is removed, 
stained, and cut into defined volumes, (d) LCM produces 
five sections per slice (outlined in yellow), by laser-melding 
onto MacroCaps (see Materials and Methods), (e) Three 
MacroCaps are used per single sample, three rats and five slices 
per Macrocap. (f) Protein eluates are produced by serially 
extracting each group of MacroCaps in the same 100 ^.L of 
extraction buffer (see text). 

(Figure le). Samples were frozen at — 70°C immo- 
bilized on the MacroCap and mounted securely into 
the opening of 0.65-mL Eppendorf tubes. A modified 
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version of the protocol, outlined by Ericsson et al. 
(2007), was used to extract whole protein from the 
immobilized brain sections. Briefly, proteins were 
extracted from the MacroCap-immobilized sections 
using sodium dodecyl sulfate-polyacrylamide gel elec- 
trophoresis (SDS-PAGE)-compatible extraction buffer 
consisting of 100 mM Tris-HCl (Sigma Aldrich, St. 
Louis, MO, USA), 0.2% Rapigest (Waters Corporation, 
Milford, MA, USA) and 1 X Protease Inhibitor Cocktail 
(Calbiochem, Billerica, MA, USA). While the samples 
were still frozen, 100 |J,L of extraction buffer per 0.13 
mm' of tissue was added to an Eppendorf tube and fit- 
ted with a MacroCap, one tube was used to pool all 
MacroCaps from each treatment group. Tubes were 
gently vortexed for 1 minute at room temperature, then 
inverted and incubated for 30 minutes at 70°C with 
shaking at 800 Hz. The tubes were then centrifuged at 
13,300 X g for 10 minutes at 4°C. The extraction from 
the first MacroCap of each treatment group was used to 
serially extract the remaining two MacroCaps, resulting 
in a pooled extract for each treatment group (Figure le). 
Samples were frozen at — 20°C and stored until further 
processing. See Supplementary Materials (available 
online at http://informahealthcare.com/doi/abs/10.3109/ 
01677063.2014.883389) for details regarding sample 
fractionation and liquid chromatography-tandem mass 
spectrometry (LC-MS/MS) analysis. 

Labeled-Pair Selection and Peptide Identification 

The LC-MS data were processed using Elucidator version 
3.2 (RosettaBiosoft, Seattle, WA) (Paweletz et al., 2010). 
Briefly, the experimental definition was built using the 
"Differential-Labeled" design type and three technical 
replicates from sample were combined and aligned using 
the PeakTeller algorithm. Labeled-pairs (correspond- 
ing heavy and light peptides) were selected using the 
Labeled-Pair Ratio Builder in Elucidator with the 
following settings: a label of K8 (8.014 Da, "C^ + ''N^), 
a maximum of three labels per peptide, 10 ppm error, and 
a retention time tolerance of ± 0.2 minutes. Raw inten- 
sity data for all "good" labeled-pairs (those meeting the 
above-mentioned criteria) were exported to Microsoft 
Excel for further analysis described below. In a sepa- 
rate analysis of the raw data in Elucidator, DTAs (text 
files corresponding to m/z and intensities from a given 
MS/MS spectra) were generated for all features and 
searched against concatenated rat forward and reverse 
sequenced database (38,401 forward protein sequences 
with entries from the Universal Protein Resource 
(www.uniprot.org). National Center for Biotechnology 
Information Reference Sequence (www.ncbi.nlm. nih. 
org), and proprietary data sets compiled as of 22 May 
2010), using MASCOT (version 2.1; Matrix Science, 
London, UK). The following search parameters were 



used: trypsin cleavage, 0.3 precursor ion and 0.8 amu 
fragment ion tolerances, variable modifications includ- 
ing methionine oxidation, acetylation of cysteine side 
chains, and K8. Peptide sequences returned from Mascot 
were further filtered at a 1% peptide false-positive rate 
using Elucidator's PeptideTeller algorithm, an extension 
of the Prophets algorithm for estimating the accuracy of 
peptide identifications (Keller et al., 2002). 

Ratio Calculation, Logic Test, and Statistical Analysis 

All ratio calculations and statistical analyses were per- 
formed using an R-language (Muenchen & Hilbe, 2010) 
script developed in-house. Labeled-pairs selected based 
on m/z and retention time as described above, and their 
respective unlabeled and labeled intensities, were exported 
from Elucidator. Labeled-pairs that returned a zero- 
intensity for the heavy-labeled peak in control samples 
(untreated labeled mixed with untreated unlabeled) were 
removed from the data set. All other zero-intensities were 
replaced with one-half of the minimum intensity (400 
units). We performed mixing error corrections, a label 
incorporation correction (as detailed in Supplementary 
Figure 1 available online at http://informahealthcare.com/ 
doi/abs/10.3109/01677063.2014.883389), and a logic test 
of the data. Briefly, unlabeled-to-labeled ratios were 
calculated for each experimental condition (observed 
ratio, R^^J and control samples (label incorporation ratio, 
R. ) based on exported intensity values. To correct for 
mixing errors, the median {R^f^J for each group was divided 
by the control R , (R ,. =0.62), resulting in 

mcorp,meaian ^ incorp,mediaii ^' ^ 

the corrected observed ratios (R ). To account for label 
incorporation differences for individual labeled-pairs, we 
divided the R by its corresponding R , resulting in 

coiT r o mcorp ^ 

the final ratio (R^^^^) as described by others (Spellman 
et al., 2008; Liao et al., 2008). The R values were then 

' ' ' norm 

subjected to a logic test. Specifically, an R value must 
be > 1.2-fold for both Ox-B and ORA groups as well 
as satisfy one of the following conditions: either (Ox-B 
ratio > Ox-B + ORA ratio > ORA ratio) or (ORA ratio 
> Ox-B + ORA ratio > Ox-B ratio). Statistical simula- 
tions were performed to estimate the false-positive rate 
of correcting and selecting labeled-pairs using the above- 
described approach. It was assumed that intensities from 
each condition have a /-distribution with a mean equal to 
that of the control (Velten, 2009). In addition, correlation 
matrices were constructed to ensure that simulated data 
(labeled and unlabeled intensities) from different con- 
ditions were comparable to the experimental data. This 
method is designed to estimate the number of peptide 
pairs falsely selected (i.e., the false-positive rate) in the 
absence of true differences between conditions (Velten, 
2009). The same procedures used for processing the 
experimental data were also applied to the simulated data 
for the calculation of R values. 
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RESULTS AND DISCUSSION 

Heavy SILAC Labeling, Weight Gain, and Litter Size 

During the observation period of 2 weeks, the four 
litters were physically inspected and weighed at days 6, 9, 
13, and 16 post weaning. All animals had the same food 
consumption and showed similar weight gain and litter 
size regardless of chow type (Figure 2). The populations 
of male littermates consisted of 15 animals with isotopi- 
cally labeled lysine and 15 control animals. To measure 
'^^Cg'^Nj-lysine incorporation into the brain proteome, 
VTA brain proteins were isolated, resolved by SDS- 
PAGE, digested in-gel, and peptides measured by high- 
resolution LC-MS. Protein labeling efficiencies were 
calculated among the labeled-pairs identified by LC-MS, 
as described in Materials and Methods and shown in 
Figure 3. The results show a distribution of heavy label 
incorporation, which is likely due to differences in rates 
of protein turnover in the brain and incomplete heavy 
label incorporation (Chevalier et al., 2009). The median 
isotopic-label incorporation efficiency among confidently 
identified proteins was 76.9%. These results are similar to 
the 71.9% of total '^N incorporation in brain in weaned 
pups by McClatchy et al. (2007). Although the incorpora- 
tion efficiency was not > 90% for all proteins, most were 
sufficiently labeled for quantitative analysis. 



Beam Break Activity Measurement 
and Compound Administration 

Phenotypic effects of the administered compounds were 
assessed by subjecting the cohorts to a beam break 
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Figure 2. Litter size and weight on normal and heavy isotope 
diet. Mean litter weights were reported as a function of the 
number of days postweaned pups were fed labeled or unlabeled 
rat chow (error bars are included but the standard deviations are 
minimal and thus not visible). HD = heavy dame (placed on the 
labeled diet 2 days into confirmed pregnancy); LD = light dame 
(fed similar, but unaltered, diet); n = the number of male pups 
included in each birth group. The data point for day 1 represents 
n = I per cohort. 
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Figure 3. In vivo heavy isotope labeling of proteins measured 
in ventral tegmental (VTA). Percent heavy label incorporation 
for labeled-pairs from '^C^'^N^-lysine-labeled Fl generation rat 
ventral tegmental tissue taken at 7 weeks post weaning. These 
data were calculated from raw MS labeled-pair intensities 
(N=297 labeled-pairs with intensity a 1E5) with median 
%Incorporation (|Ij,2) = 76.9%. Heavy label incorporation 
calculations predicted > 100% incorporation for 22 labeled- 
pairs (average: 108%), demonstrating a small contribution of 
error to labeled-pair ratio calculations likely from mixing and 
changes in isotopic distributions between heavy and light amino 
acid-containing peptides. 



behavioral data assay. The data show the orexin receptor 
antagonist was effective at lowering basal locomotor 
activity relative to vehicle prior to Ox-B administration 
(Figure 4). As previously demonstrated, animals treated 
with Ox-B showed a marked increase in activity relative 
to control (Winrow et al., 2010). Pretreatment with the 
orexin receptor antagonist significantly inhibited the 
Ox-B-induced locomotor response (Figure 4). Animals 
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Figure 4. Orexin receptor antagonist response based on 
beam break behavioral data. Locomotor activity in response to 
orexin (Ox-B) and orexin receptor antagonist (ORA). Animals 
were treated with control (vehicle) or ORA at 10:30 am 
(as indicated, A) and treated with aCSF or orexin peptide 
at 11:30 AM (as indicated, T). Significant ORA differences 
between control (vehicle) vs. Ox-B and ORA vs. Ox-B 
(*P value < 0.05) and between ORA vs. Ox-B and ORA vs. 
aCSF (*P value < 0.05) are shown. 
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administered Ox-B an hour after ORA administration 
showed a modulated response, resulting in -500 beam 
breaks above baseline of Ox-B administration (-1800 
beam breaks), demonstrating that the antagonist effec- 
tively inhibited Ox-B-induced locomotor activity. 



Laser Capture Microdissection of Proteins 
From the VTA, Labeled-Pair Selection, 
and Peptide Identification 

The primary purpose of performing differential pro- 
teomic analysis on a very specific subregion of the brain 
in these studies is to reduce the tissue heterogeneity and 
to increase the likelihood of identifying proteins related 
to sleep biology. The disadvantage of this approach is 
that a more complicated sample processing is required 
and protein recovery is greatly reduced if compared with 
whole brain protein homogenates. The Elucidator results 
returned 44,163 detected features (with specific retention 
times, ni/z, and intensity values), which collapsed into 
28,354 isotope groups. Labeled-pairs were selected in a 
"counterclockwise" fashion whereby the specific mea- 
surable characteristics (m/z and retention time, criteria 
described in Materials and Methods) of individual isotope 
groups were used to establish links between heavy- and 
light-labeled species, resulting in the selection of 239 
labeled-pairs. This approach is deemed "counter" to the 
typical "clockwise" approach to SILAC data where pairs 
are determined based on peptide identification, making 
confidently assigned tandem mass spectra a prerequisite 
for any analysis of labeled-pairs. The counterclockwise 
approach allows for selection of pairs based on accurate 
mass with the quantitative characteristics of interest, inde- 
pendent of peptide identification. A separate data analysis 
was performed in Elucidator to determine what proteins 
were identified via data-dependent MS/MS (shot-gun) 
in this experiment, resulting in the identification of 
548 proteins. Additional targeted analyses were performed 
to generate tandem mass spectra (MS/MS) for label pairs 
of interest that were not generated by shot-gun MS/MS 
analysis. 

In total, as shown in Table 1, 15 labeled-pairs 
were determined to be of interest either (Ox-B ratio 

> Ox-B + ORA ratio > ORA ratio) or (ORA ratio 

> x-B + ORA ratio > Ox-B ratio) based on the data treat- 
ment and criteria described in Materials and Methods 
and Supplementary Figure 1 available online at http:// 
informahealthcare.com/doi/abs/1 0. 3 1 09/0 1 677063 .20 14. 
883389 (false-positive rate estimated to range from 16% 
to 38% with n = 500 simulations and 2-20 degrees of free- 
dom as described in Materials and Methods). Eleven of 
the 15 selected pairs had tandem mass spectra present in 
the original experiment that allowed for confident peptide 
identification, 2 were sequenced with a targeted approach 
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using inclusion lists of m/z values, and 2 remain unidenti- 
fied due to insufficient quality MS/MS spectra. Among the 
13 identified pairs of interest, 7 increased in relative abun- 
dance in response to orexin receptor antagonist treatment 
and were identified as peptides derived from; L-lactate 
dehydrogenase A (LDHA), Na/K-transporting ATPase 
(ATPIBI), 2-,3'-cyclic nucleotide 3 '-phosphodiesterase 
(CNP), synaptotagmin (SYTl), tubulin (TUB), and 
14-3-3ri (YWHAH). Six labeled-pairs decreased in rela- 
tive abundance as a result of antagonist treatment and were 
identified as peptides from myelin basic protein (MBP), 
mitochondrial malate dehydrogenase, peptidyl-prolyl 
cisltrans isomerase A (PPIA), a-synuclein (SNCA), and 
ADP-ribosylation factor 3 (ARFl). Notably, two proteins 
were identified by two individual peptides each (Table 1). 
Both labeled-pair ratios for PI 3233 (CNP) showed down- 
regulation upon ORA treatment, with higher variability 
with respect to Ox-B treatment (coefficient of variabil- 
ity [CV] = 81%) than Ox-B ORA (20.4%) and ORA 
(33.6%); however, all ratios trended in the same direction. 
PlOlll (PPIA) was up-regulated similarly with regard 
to Ox-B treatment (CV = 0.6%), and closely tracked in 
both Ox-B+ORA (3.5%) and ORA (15.7%) treatments. 
Supplementary Figures 2 and 3 available online at http:// 
informahealthcare.com/doi/abs/10.3109/01677063. 
2014.883389 contain the MS/MS spectra for peptides 
presented in Table 1 . 

Recent findings by Waanders et al. (2009) showed 
that as few as 2000-4000 cells can provide a sufficient 
number of cells for protein differential analysis in pan- 
creatic islets. In addition, recently published LCM- 
LC-MS/MS data, using a similar workflow to the one 
presented here, showed that it is feasible to identify and 
differentially quantitate a significant number of proteins 
from microdissected neuronal tissue (Han et al., 2008; 
Waanders et al., 2009). Due to regional volume con- 
straints in the VTA in these studies, two 13-mm^ tissue 
sections (average area as calculated by the Veritas LCM 
software; Molecular Devices) were captured per animal 
(Figure 2c-e), with three animals pooled per condi- 
tion. Based on a total tissue volume of 0.39 mm' with 
a mean eukaryotic cell diameter of 20 jim, we estimate 
the cell count at approximately 180,000 per condition. 
It is likely that a small portion of this disparity between 
the number of proteins identified in these experiments 
versus the number identified by Waanders et al. (2009) 
can be attributed to obvious cell type difference between 
neurons and islets, primarily the large membrane to cell 
volume ratio in neurons resulting from highly com- 
plex cellular process structure as well as myelination. 
In addition, in this procedurally complex, multiarm 
study several important factors seem to be at the cen- 
ter of ensuring overall reproducibility and reducing 
sample loss during biochemical processing. First, in 
vivo isotopic labeling of Fl generation rats results in 



incomplete label incorporation due to different turnover 
rates found among the protein population. Continuous 
labeling into a second F2 generation cohort have dis- 
tinct benefits as reported by McClatchy et al. (2007), 
resulting in increased label incorporation in neuronal 
tissue. These studies, however, were limited by the time 
lines associated with labeling a second generation of 
rats. Second, consistent VTA sample volumes were not 
strictly obtainable (usually ± 1 mm'). Specific normal- 
ization measures were incorporated to correct for this 
error by standardizing protein extraction buffer volume 
as a function of reported tissue volume; however, pre- 
cise cell populations may not have been reproducibly 
solubilized, thereby perturbing intercellular peptide and 
protein concentrations between samples. The protein 
extraction buffer used in this experiment was not condu- 
cive to conventional protein concentration assays, and 
thus concentrations were normalized by the tissue 
volume provided by LCM calculations. Third, the lysis of 
neuronal tissue sections physically bound to MacroCap 
surfaces may be inconsistent between samples due to 
the lack of physical agitation (e.g., sonication or hyper- 
tonic disruption). Therefore, certain proteins, especially 
those bound to membranes, may not have been con- 
sistently represented. The extraction buffer used was 
optimized for total protein content and "distribution" 
along the lane as revealed by Coomassie blue stain- 
ing on an SDS-PAGE gel (data not shown), but it too 
could be further optimized or expanded to perform 
stepwise solubility improvement and/or cellular fraction 
concentration. 



Biological Relevance of Protein Identifications 

To understand the functional significance of the altera- 
tions in protein expression. Ingenuity Pathway Analysis 
(Ingenuity Systems, Redwood City, CA) was used to 
interrogate the identified proteins and to identify over- 
lapping associations with known orexin signaling. These 
observed associations, though intriguing, will benefit 
from further independent validation and are affected 
by the limited number of peptides identified in the cur- 
rent study. This evaluation revealed that key proteins 
involved in mammalian neuron growth, development, 
and energy metabolism were altered in response to the 
treatment conditions in the VTA along with several 
proteins with defined roles in synaptic function, behav- 
ior, and locomotion (Figure 5). Additionally, many of 
these proteins have been associated with neurological 
disorders (ATPIBI, CNR LDHA, MBP, PPIA, SNCA, 
SYTl, TUBAIA, YWHAH) (Supplementary Table 1 
available online at http;//informahealthcare.com/doi/ 
abs/10.3109/01677063.2014.883389). A number of 
associations with neurodegenerative and psychiatric 
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Figure 5. Ingenuity network analysis of differentially regulated labeled-pairs, proteins listed in Table 1 were evaluated using Ingenuity 
Pathway Analysis (version 7.6) to determine interacting nodes. The top network identified was "Neurological Disease, Cell Death, 
Nervous System Development and Function" and is shown along with the L/H normalized ratios for each protein. Proteins with 
significant L/H ratios for Ox-B treatment are indicated in red and those with significant L/H ratios for antagonist treatment are shown 
in green, with ratios values shown below. 



disorders were identified, including schizophrenia 
(CNP, SYTl, TUBIA, YWAH), Parkinson's disease 
(LDHA, MBP, SNCA, TUBAIA), Huntington's dis- 
ease (LDHA, PPIA, TUBAIA), Alzheimer's disease 
(SNCA), and multiple sclerosis (MBP). Most of the 
proteins listed (10/11) are localized in the cytoplasm, 
with functions including synaptic transmission (ARFl, 
CNP MBP SNCA, SYTl), axonal morphology (CNP 
MBP, SNCA, YWHAH), and neurotransmitter signaling 
(CNP SNCA, SYTl). Importantly, many of the differ- 
entially expressed proteins have intersecting functional 
associations with orexin signaling. These include roles 
in synaptic transmission, neurotransmitter release, and 
neuronal survival. Further, given the clear increase in 
beam break activity following administration of Ox-B 
peptide, it is interesting to observe that CNP, MBP, and 



SNCA levels were altered, as these proteins are involved 
in neuronal activation and locomotor activity. 

Although there were no direct physical connections 
identified between orexin receptors, peptides, and the 
altered proteins, the Ingenuity nodal analysis revealed 
associated overlapping ontological pathways, including 
those involved with the clock-controlled cell cycle gene 
Myc, which has shown circadian periodicity in cancer 
studies (lurisci et al, 2006; Levi et al., 2007). MBP, active 
in the myelination of nerves, is also associated with Myc 
functioning outside of the cellular space. MBP and SNCA 
(the "non-AP" component of Alzheimer's disease) are 
both down-regulated by the antagonist affecting a 
number of pathways, including central nervous system 
(CNS) morphology, nervous tissue cell viability, CNS 
development, and fibril formation. In addition, MBP and 
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SNCA appear to have an opposing expression relation- 
ship with CNP (ambiguously associated with myelin) 
and YWHAH (involved in various signaling pathways) in 
neurite development and human neurological disorders. 

CNP is also oppositely modulated compared with MBP 
under the influence of antagonist treatment. The proteins 
are conversely expressed during neuron and axon develop- 
ment, filament assembly, and, as appropriate for this study, 
locomotion. The alpha-l-A chain of tubulin (TUBAIA) 
interacts with both CNP and MBP (Figure 5) and is up- 
regulated with regard to receptor antagonist treatment. The 
functions of SNCA continue to emerge, and the question 
remains as to why CNP remains active during induced sleep 
whereas MBP appears to be decreased by the ORA treat- 
ment. The change in SNCA levels in response to Ox-B in 
the VTA is intriguing given the recent linkage of this gene 
with addictive behaviors in rodents (Spence et al., 2009) 
and humans (Clarimon et al., 2007; Bonsch et al., 2005). 
Although many of the peptides identified have been linked 
with CNS function and some pathologies, the relevance 
of the observed changes in VTA in response to peptide or 
antagonist treatment remains to be determined. 

SUMMARY 

Several recent studies have shown the key role that orexin 
signaling plays in addictive behaviors modulated through 
the VTA, including the ability to block amphetamine sen- 
sitization and nicotine reinstatement with selective orexin 
receptor antagonists (Winrow et al., 2010; Aston-Jones 
et al., 2009, 2010). The studies described here strengthen 
the role for orexin signaling through the VTA and point 
to a previously unknown connection between SNCA and 
orexin signaling in VTA function. It will be important to 
further confirm the peptide changes observed in the cur- 
rent study, and additional investigations will be required 
to understand the potential relevance of these changes 
in rodent VTA. In addition, improving the proteomics 
methodology to target an increased and a more diverse 
total protein extraction from LCM-isolated tissue is a 
reasonable expectation to strive for in future studies. 
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